1/1 


TO-H140  439 
UNCLASSIFIED 


TRANSPORT  OF  HATER  DISSOLVED  SUBSTANCES  HEAT  AND 
ELECTRIC  CURRENT  THROUGH.  .  (U>  ILLINOIS  UNIV  AT  URBANA 
DEPT  OF  GEOLOGV  D  L  GRAF  ET  AL.  20  HAR  84 
RRO-16743.  S-GS  OAAG29-80-C-0068  F/G  8/13 


NL 


MICROCOPY  RESOLUTION  TEST  CHART  f 
NATIONAL  BUREAU  OF  STANDARDS-1963-A 


/ffio  /^7¥3.S 


TRANSPORT  OF  WATER,  DISSOLVED 
SUBSTANCES,  HEAT  AND  ELECTRIC 
CURRENT  THROUGH  SHALES  AND 
CLAY-RICH  SEDIMENTS 


FINAL  REPORT 


DONALD  L.  GRAF 
WILLIAM  M.  BENZEL 
ILHAM  DEMIR 
PAUL  HAYDON 


MARCH  20,  1984 


U.  S.  ARMY  RESEARCH  OFFICE 


CONTRACT  DAAG  29-80-C-0068 


DEPARTMENT  OF  GEOLOGY 
UNIVERSITY  OF  ILLINOIS 
AT  URBANA  -  CHAMPAIGN 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


W-, 


'.•*  V  V  * 


Unclass i f ied 


1  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

I.  NEPORT  number 

ARO  16743.5-GS 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  catalog  number 

4.  Title  (and  SubfIlU) 

Transport  of  Water,  Dissolved  Substances,  Heat 
and  Electric  Current  Through  Shales  and 

Clay-Rich  Sediments 

5.  type  of  report  a  period  covered 

1  Feb  1980  -  31  Jan  1984 

Final  Report 

s.  performing  org.  report  number 

T.  AUTHORfa; 

Donald  L.  Graf  Paul  llaydon 

William  M.  Benzel 

Ilham  Demir 

B.  contract  or  grant  numberc*) 

DAAG29-80-C-0068 

>■  PERFORMING  ORGANIZATION  NAME  ANO  AOOnrSS 

Illinois,  Univ 

10.  program  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 

N/A 

H.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

U.  S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 

12.  report  date 

13.  number  OF  PAGES 

13 

U.  MONITORING  AGENCY  NAME  A  ADDRESSfl/  dlltottnl  Iram  ControUIng  Ollle») 

IS.  SECURITY  CLASS,  fof  thie  report) 

Unclassified 

ISa.  declassification/ DOWNGRADING 
SCHEDULE 

IS  distribution  STATEMENT  rot  thit  Rtpotl) 

Approved  for  public  release;  distribution  unlimited. 

IT.  DISTRIBUTION  STATEMENT  (al  lh»  abalraef  mtrrr*d  In  Black  70.  It  dlllcrcnt  hmn  Report) 

IS.  supplementary  notes 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  , Policy,  or  decision,  inless  so  designated  by  other  documentation 

19.  KEY  WORDS  ('Conflnuw  on  f«r«r«o  o/do  it  noco«««ry  mnd  id^nUty  by  6locJ^  numb9r) 

Water  Transport  Properties  Electric  Current 

Sediment  Heat  Transfer  Membranes 

Clay  Temperature 

Research  effort  was  concentrated  upon  building  two  heavy-duty 
membrane  presses,  securing  adcpiate  bulk  samples  of  natural  clays,  and  reviewing 
literature.  The  principal  portion  of  the  contract  effort  reported  here  is  con¬ 
tained  in  two  completed  Ph.D.  theses,  by  William  M.  Benzel  and  Paul  R.  Haydon, 


rOMM 


M)  I  jam  TS 


EOlTtOM  OF  I  MOV  6S  IS  OUSOLETC 


- 1473  - - UNCLASSIFIED _ 


Unclassified 


security  classification  of  this  PACE(Whan  Dai«  BnUnd) 


ARO  16743.5-GS 
20.  ABSTRACT  CONTINUED: 


and  one  that  is  nearing  completi'ti,  by  Ilham  Demlr.  These  theses  consider  the 
effect  upon  clay-membrane  transpi'rt  at  20“  C  of  the  thickness  and  fabric  of  the 
clay  cake,  the  effect  of  temperature  in  the  range  20“-180“  C  upon  membrane 
transport  and  membrane  fractionation  of  the  stable  oxygen  isotopes  contained  in 
water  molecules,  and  the  electmkinetic  behavior  of  these  clay  membranes  at  20“ 
C.  It  becomes  possible  after  c<Miipletlon  of  these  many  months  of  experimental 
runs  to  evaluate  the  difficultli  ;  encountered  with  experimental  apparatus. 


security  classification  of  this  PACE(T»h»fi  Data  Bnlararf) 


TRANSPORT  OF  WATER.  DISSOLVED 
SUBSTANCES.  HEAT  AND  ELECTRIC 
CURRENT  THROUGH  SHALES  AND 
CLAY-RICH  SEDIMENTS 


FINAL  REPORT 


DONALD  L.  GRAF 
WILLIAM  M.  BENZEL 
ILHAM  DEMIR 
PAUL  HAYDON 


MARCH  20.  198A 


U.  S.  ARMY  RESEARCH  OFFICE 


CONTRACT  DAAG  29-80-C-0068 


1  j  on  For 


DEPARTMENT  OF  GEOLOGY 
UNIVERSITY  OF  ILLINOIS 
AT  URBANA  -  CHAMPAIGN 


Vx  ) 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


-.'o 


THE  FINDINGS  IN  THIS  REPORT  ARE 
NOT  TO  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  ARMY  POSITION. 
UNLESS  SO  DESIGNATED  BY 
OTHER  AUTHORIZED  DOCUMENTS. 


TABLE  OF  CONTENTS 


INTRODUCTION 

EFFECT  UPON  TRANSPORT  OF  VARYING 
THICKNESS  AND  FABRIC  OF  CLAY 
LAYER 

TEMPERATURE  DEPENDENCE  OF  TRANSPORT 

FRACTIONATION  OF  STABLE  ISOTOPES 
DURING  TRANSPORT 

ELECTRICAL  ASPECTS  OF  TRANSPORT 

RELIABILITY  AND  PERFORMANCE  OF 
EXPERIMENTAL  APPARATUS 

A  MODEL  FOR  AQUEOUS  ELECTROLYTE 
DIFFUSION 

IMPORTANCE  OF  REVERSE  CHEMICAL 
OSMOSIS  IN  NATURE 

MINERAL  DEHYDRATION  IN  SEDIMENTARY 
ENVIRONMENTS 

PUBLICATIONS  AND  ORAL  PRESENTATIONS 


INTRODUCTION 

Iff 

Contract  DAAG  29-80-C-0068  is  a  continuation  of  Grant  DAAG^77-G-0011 . 
Research  effort  during  the  grant  was  concentrated  upon  building  two  heavy-duty 
membrane  presses,  securing  adequate  bulk  samples  of  natural  clays,  and  reviewing 
literature.  The  principal  portion  of  the  contract  effort  reported  here  is  con¬ 
tained  in  two  completed  Ph.D.  theses,  by  William  M.  Benzel  and  Paul  R.  Haydon, 
and  one  that  is  nearing  completion,  by  Ilham  Demlr.  These  theses  consider  the 
effect  upon  clay-membrane  transport  at  20**  C  of  the  thickness  and  fabric  of  the 
clay  cake,  the  effect  of  temperature  in  the  range  20**-180‘*  C  upon  membrane 
transport  and  membrane  fractionation  of  the  stable  oxygen  Isotopes  contained  in 
water  molecules,  and  the  electroklnetic  behavior  of  these  clay  membranes  at  20** 
C.  It  becomes  possible  after  completion  of  these  many  months  of  experimental 
runs  to  evaluate  the  difficulties  encountered  with  experimental  apparatus. 

The  detailed  mechanisms  involved  in  transport  through  the  micropores  of  the 
compacted  clay  cake  are  of  considerable  scientific  interest.  Donald  L.  Graf 
spent  some  contract  project  time  developing  a  model  for  electrolyte  diffusion  in 
bulk  aqueous  solutions,  as  a  preliminary  to  understanding  the  more  complex  dif¬ 
fusion  processes  within  the  pores.  Mindful  of  the  skepticism  among  many  geolo¬ 
gists  about  the  significance  of  reverse  chemical  osmosis  in  nature,  he  evaluated 
the  existence  in  nature  of  fluid  pressure  differentials  adequate  to  drive  this 
process.  Together  with  David  E.  Anderson,  he  reviewed  the  pressure,  tem¬ 
perature,  and  salinity  dependences  of  mineral  dehydration  reactions  that  could 
release  liquid  water  in  thick  sedimentary  sections. 

These  research  results  are  presented  in  detail  in  theses  and  in  scientific 
papers  that  have  either  been  published  or  are  in  preparation.  This  report 
therefore  limits  itself  to  offering  qualitative  summaries  of  the  principal  fin¬ 
dings. 


EFFECT  UPON  TRANSPORT  OF 
VARYING  THICKNESS  AND  FABRIC 
OF  CLAY  LAYER 

These  experiments  were  run  at  20°  C  under  a  confining  pressure  of  5000  psi, 
using  a  solution  0.79  molal  in  NaCl  and  0.079  molal  in  CaCl^*  The  relatively 
low  salinity,  relative  to  the  total  range  observed  in  nature,  was  chosen  for 
ease  of  comparison  of  results  with  the  published  engineering  literature  on  near¬ 
surface  slope  stability  problems. 

The  clay  cakes  used  were  constructed  by  dispersing  the  clay  crystallites  in 
water,  separating  out  a  narrow  range  of  crystallite  sizes  (0.5  to  2.0  ym 
diameter),  and  then  sedimenting  these  particles.  The  labor  Involved  made  it 
impractical  to  construct  really  thick  cakes.  Within  the  fourfold  range  from  0.5 
to  2.0  cm,  thickness  appears  to  have  little  or  no  effect  upon  filtration 
efficiency. 

It  was  possible  by  changing  the  salinity  of  the  solution  to  form  a  cake  with 
either  relatively  random  or  relatively  well-oriented  fabric.  These  fabrics 
could  be  substantially  preserved  if  compaction  was  carried  out  slowly  and  care¬ 
fully.  The  salt-filtration  efficiencies  of  the  cakes  with  oriented  fabric  are 

roughly  twice  as  great  as  those  of  the  others,  reaching  50  to  68%.  There  was 
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preferentlal  passage  of  Ca  through  the  cake,  relative  to  Na  ,  in  three  of  the 


f  our  runs 


TEMPERATURE  DEPENDENCE  OF  TRANSPORT 

These  experiments  were  run  at  20*,  95®,  140®  and  180®  C  under  confining 
pressures  of  5000  and  10,000  psi,  using  clay  cakes  with  random  fabric  and  a 
solution  5.0  molal  in  NaCl  and  0.45  molal  in  CaCl^*  The  temperature  range  and 
the  high  salinity  were  Intended  to  simulate  depths  of  5000  to  20,000  feet  in 
thick  sedimentary  sections. 

The  filtration  efficiencies  for  these  experiments  range  from  16  to  28%  for 
NaCl  and  3  to  32%  for  CaCl^ •  Those  for  CaCl2  Increase  nearly  linearly  with  tern 
perature,  whereas  those  for  NaCl  go  through  a  minimum  at  the  two  intermediate 
temperatures.  The  crossover  temperature,  above  which  CaCl^  is  filtered  more 
efficiently  than  NaCl,  is  about  160®  C  at  5000  psi  compaction  pressure.  At 
10,000  psi,  it  is  somewhere  below  140®  C.  , 

Qualitative  arguments  involving  the  pressure  and  temperature  dependences  of 
various  physical  parameters  suggest  that  fluid  viscosity  must  be  involved  in 
explaining  the  general  Increase  in  filtration  efficiency  with  temperature.  The 
effect  of  compaction  pressure  upon  mean  pore  size  affords  a  satisfactory  expla¬ 
nation  for  the  Increased  filtration  efficiencies  observed  at  the  higher  compac¬ 
tion  pressure. 

Quantitative  calculations  using  the  Teorell-Meyer-Slever  and  Fritz-Marine 
models  for  predicting  effluent  chemical  composition  are  partially  completed. 

The  TSM  values  at  20®  C  are  in  excellent  agreement  with  experimental  values . 


FRACTIONATION  OF  STABLE  OXYGEN  ISOTOPES 

The  isotopic  enrichment  of  the  oxygen  in  water  passing  through  a  well- 
oriented  smectite  cake  at  20°  C  is  0.96  o/oo,  in  good  agreement  with  the  value 
of  -0.8  o/oo  found  earlier  by  T.  B.  Coplen  and  B.  B.  Hanshaw.  Possible  small 
shifts  in  this  value  occasioned  by  change  of  fabric  or  cake  thickness  are  masked 
by  uncertainties  In  our  experimental  procedures. 

The  enrichments  found  at  140“  and  180“  C,  0.15  and  0.05  o/oo,  are  effec¬ 
tively  zero  considering  the  uncertainties  In  measurement  and  data  reduction.  We 
are  attempting  to  decide  whether  the  rate  of  decrease  of  enrichment  between  20“ 
and  140“  C  favors  a  particular  one  of  the  possible  Isotopic  fractionation  mecha¬ 


nisms 


ELECTRICAL  ASPECTS  OF  TRANSPORT 

These  experiments  are  being  run  at  20°  C  under  a  confining  pressure  of  5000 
psl,  using  clay  cakes  with  either  random  or  oriented  fabric  and  either  a  1.087  m 
NaCl  solution  or  a  solution  0.92  m  In  NaCl  and  0.075  m  In  CaCl^*  The  membrane 
press  Is  modified  for  these  electroklnetlc  studies  by  adding  four  planar  Ag-AgCl 
mesh  electrodes,  two  on  either  side  of  the  clay  cake.  In  one  type  of 
experiment,  for  example,  a  voltage  potential  Is  Imposed  across  a  pair  of 
electrodes  (one  on  either  side  of  the  cake)  and  the  flow  of  current  Is  monitored 
using  the  other  pair. 

Streaming  potential  values  decreased  by  18  to  27%  during  the  8  to  10  weeks 
It  took  for  mechanical  and  chemical  steady  state  to  be  reached.  This  deter¬ 
mination  affords  sode  Insight  Into  the  effect  that  continuing  compaction  In 
nature  might  have  on  values  obtained  In  borehole  logging.  Varying  hydraulic 
potential  produced  an  approximately  linear  variation  In  streaming  potential,  but 
the  plot  of  these  two  variables  had  a  different  slope  for  each  run. 

Hydraulic  conductivities  found  In  runs  In  which  the  streaming  potential  Is 
allowed  to  develop  are  about  one-fifth  those  for  runs  In  which  streaming  poten¬ 
tial  Is  short-circuited,  so  that  streaming  potential  and  the  electrovlscous 
effect  modify  significantly  the  hydraulic  transport  through  shales. 

Adding  CaCl^  In  modest  amounts  to  the  NaCl  solution  has  an  appreciable 
effect.  The  Ca  dominates  electroklnetlcally  over  the  Na.  It  appears  that  the 
potential  difference  at  the  downstream  cake-brine  Interface  Is  mainly  respon¬ 
sible  for  salt  filtration  In  the  NaCl  system,  but  that  the  upstream  Interface  Is 
equally  Important  In  the  NaCl-CaCl2  system. 

Saxen's  Law  of  nonequilibrium  thermodynamics,  which  relates  the  magnitudes 
of  current  and  fluid  flows  for  Impose  fluid  and  electrical  potentials,  Is 


obeyed  very  closely  by  the  experimental  system  smectlte-1 .087  m  NaCl  solution  at 
room  temperature.  This  conclusion  was  made  possible  by  a  new  data  analysis  pro¬ 
cedure  that  corrects  for  such  effects  as  fluid  compressibility  and  chemical 
osmotic  transport. 


RELIABILITY  AND  PERFORMANCE  OF 
EXPERIMENTAL  APPARATUS 


The  membrane  presses  have  required  a  good  deal  of  machine-shop  and 
electronic-shop  support  for  repair  and  minor  modification.  The  amount  of  down 
time  has  not  been  excessive  for  such  complex  apparatus,  but  It  has  generally  not 
been  possible  to  sustain  experiments  of  several  months  without  Interruption. 
Valving  In  the  system  has  made  It  possible  to  preserve  the  clay  cake  and  Its 
Immediately  surrounding  fluid  environment  during  most  Interruptions,  but  pertur¬ 
bations  of  fluid  composition  In  the  Input  and  stirring  reservoirs  have  extended 
the  time  to  reach  steady  state. 

Specific  difficulties  have  Included  leaks  because  of  scoring  of  moving 
pistons,  requiring  polishing  of  these  components,  cracking  of  the  hlgh-purlty 
alumina  liner  In  the  sample  assembly  when  temperature  was  raised  high  enough  to 
release  the  compressive  stress  placed  upon  the  liner  earlier  by  shrink  fitting, 
and  deterioration  of  Vlton  and  ethylene  polypropylene  0-rlngs  because  of  heat, 
abrasion,  and  possibly  chemical  attack.  The  advantages  of  Teflon  0-rlngs  In 
resisting  temperature  and  chemical  attack  could  not  be  utilized  because  this 
material  Is  so  stiff  that  It  did  not  deform  In  response  to  changing  pressure, 
and  leaks  resulted. 

A  certain  amount  of  leakage  can  be  tolerated,  by  Increasing  the  fluid  Input 
rate  to  compensate.  Mass  balance  calculations  are  made  more  difficult,  but 
they  are  not  very  satisfactory  anyway  because  of  uncertainty  about  how  much 
water  Is  expelled  from  the  clay  cake  during  the  final  stages  of  compaction. 


A  MODEL  FOR  AQUEOUS  ELECTROLYTE  DIFFUSION 


Current  chemical  models  for  mutual  diffusion  in  aqueous  electrolyte  solu¬ 
tions  derive  from  the  classic  work  of  Lars  Onsager  in  the  1930' s.  Although  his 
treatments  of  the  electrophoretic  and  relaxation  effects  are  Impressive  mathema 
tlcal  feats,  the  failure  of  Onsager-type  models  to  predict  diffusion  coef¬ 
ficients  at  concentrations  beyond  about  0.05  molar  suggests  that  the  physical 
assumptions  Involved  are  Inadequate  or  erroneous.  The  model  proposed  by  D.  L. 
Graf,  D.  E.  Anderson  and  J.  B.  Uoodhouse  starts  by  unifying  currently  available 
descriptions  of  diffusion  based  upon  nonequilibrium  thermodynamics  and 
transition-state  theory.  The  residual  error  between  the  predictions  of  this 
combined  theory  and  the  values  obtained  by  experimental  measurement  can  be 
largely  eliminated  by  Introducing  a  negative  exponential  In  concentration. 
Adjustment  of  the  single  variable  parameter  In  the  exponential  gives  diffusion- 
coefficient  prediction  for  NaCl-H20,  KCI-H2O  and  Na2SO^-H20  within  the  experi¬ 
mental  error  of  measurement,  0.1-0. 2%.  The  maximum  error  of  prediction  for 
Individual  compositions  within  CaCl2“H20,  SrCl2~H20,  BaCl2“H20,  MgCl2“H2® 
MgS0^-U20  Is  about  4%,  a  huge  Improvement  over  what  was  possible  previously. 

The  experimentally-measured  diffusion  coefficients  can  be  reproduced  with 
errors  comparable  to  those  of  the  empirical  fits  by  further  postulating  that 
Individual  Ion-water  molecule  exchanges  are  coupled  to  yield  hydrated  neutral 
exchange  complexes  (the  activated  complexes),  and  producing  the  required  nega¬ 
tive  expotentials  by  calculating  probabilities  for  specific  complexes. 
Verification  by  other  experimental  methods  of  at  least  one  of  the  hydrated 
neutral  exchange  complexes  would  be  desirable,  but  It  appears  that  not  a  great 
deal  Is  known  about  the  hydration  state  of  such  neutral  species. 


IMPORTANCE  OF  REVERSE  CHEMICAL 
OSMOSIS  IN  NATURE 


A  number  of  workers  In  the  decades  following  World  War  II  have  Invoked 
reverse  chemical  osmosis  to  help  explain  subsurface  brine  compositions.  Water 
moves  upward  across  shales  from  the  more  saline  brines  below,  overcoming  the 
downward-directed  chemical  osmotic  pressure  differential.  Doubts  about  this 
mechanism  Involve  the  existence  In  nature  of  adequately  large  upward-directed 
fluid  pressure  gradients,  and  the  possibility  that  the  observed  brine  com¬ 
positions  might  more  easily  be  explained  by  Ion  exchange  and/or  chemical  reac¬ 
tion. 

Donald  L.  Graf  concluded  that  an  unequivocal  demonstration  of  the  operation 

of  reverse  osmosis  In  nature  probably  requires  that  several  parameters,  e.g., 

4*  2*4 

Na  concentration,  Ca  concentration,  fluid  pressure,  and  self-potential  be 
shown  to  have  the  proper  values  to  have  resulted  from  the  effects  of  a  shale 
membrane  superposed  on  the  regional  values  of  these  quantities. 

Calculations  showed  that  overpressuring  from  the  rapid  deposition  of  fine¬ 
grained  sediments  was  more  than  adequate  to  drive  reverse  chemical  osmosis,  but 
that  from  topographic  relief  was  not.  Occurrence  of  overpressuring  In  sedimen¬ 
tary  sections  older  than  Cretaceous  indicate  that  post-depositional  mechanisms 
such  as  tectonic  compression  and  aquathermal  pressuring  must  also  operate. 

In  published  discussion  after  the  appearance  of  Graf's  paper.  It  was  pointed 
out  that  he  had  not  considered  the  Inefficiency  of  chemical  osmosis,  and  that 
even  topographic  relief  was  adequate  in  many  localities  to  drive  reverse  chemi¬ 
cal  osmosis. 


9 


MINERAL  DEHYDRATIONS  IN  SEDIMENTARY  ENVIRONMENTS 

Calculations  about  the  effect  of  shale  membrane  processes  upon  fluid 
pressure  and  brine  composition  distributions  In  subsurface  sedimentary  environ¬ 


ments  have  to  consider  the  water  added  to  the  system  with  time  because  of 
mineral  dehydrations. 

D.  L.  Graf  and  D.  E.  Anderson  reviewed  what  was  known  about  the  dehydrations 
of  gypsum,  smectite,  halloyslte,  vermlcullte,  and  the  zeolite  minerals.  Simple 
dehydrations  such  as  those  of  gypsum  and  halloyslte  occur  at  sharply-defined 
temperatures  and  thus  contribute  a  tlme-llmlted  fluid  pulse  at  a  given  point. 

The  dehydrations  of  analclme  and  smectite  proceed  by  reactions  Involving  other 
sedimentary  minerals.  The  smectite  reaction  proceeds  by  way  of  a  succession  of 
mixed-layered  llllte/smectltes  and  may  require  several  million  years  for  comple¬ 
tion  at  a  given  stratigraphic  Interval.  The  great  variety  of  zeolite  com¬ 
positions  and  crystal  structures  makes  It  difficult  to  generalize  about  their 
behavior . 
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PUBLISHED  PAPERS  AND  ORAL  PRESENTATIONS 


Because  most  of  D.  L.  Graf's  files  are  still  in  storage  after  moving  from 
Illinois  to  Arizona,  it  has  been  necessary  to  give  paraphrased  titles  for  some 
of  the  talks  listed  belov.  These  titles  are  shown  In  parentheses. 

Oral  Presentations 

Graf,  D.  L.,  and  Anderson,  D.  E.,  (Contributions  from  diffusion  and  membrane 
transport  to  subsurface  water  compositions):  Invited  paper.  Joint  Meeting  of 
the  International  Geological  Congress  and  the  International  Association  of 
Scientific  Hydrology,  Paris,  July,  1980. 

Benzel,  H.  M. ,  and  Graf,  D.  L.,  (Experimental  measurement  of  transport 
through  clay  layers  at  room  temperature):  Clay  Minerals  Society  Meeting, 
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